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1 Abstract 

We present self-consistent calculations for the self-energy and magnetic sus- 
i ■ ceptibility of the 2D and 3D symmetric Anderson lattice Hamiltonian, in the 

g 

o 

o 



fluctuation exchange approximation. At high temperatures, strong f-electron 
scattering leads to broad quasiparticle spectral functions, a reduced quasi- 



^ ■ particle band gap, and a metallic density of states. As the temperature is 

lowered, the spectral functions narrow and a pseudogap forms at the char- 
acteristic temperature T x at which the width of the quasiparticle spectral 
function at the gap edge is comparable to the renormalized activation energy. 
For T << T x , the pseudogap is approximately equal to the hybridization gap 
in the bare band structure. The opening of the pseudogap is clearly apparent 
in both the spin susceptibility and the compressibility. 
PACS numbers: 71.28.+d, 75.20.hr 
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It is well known that heavy electron systems may be superconducting, magnetic or para- 
magnetic at the lowest temperatures Recent discoveries of cerium and uranium based 
compounds with quasiparticle gaps or pseudogaps E g ~ 10 — 100K have generally been 
interpreted as evidence that the zoo of heavy electron systems includes semiconductors as 
well . However, the interpretation of the data is not unambiguous, and it is not yet clear 
whether these systems should be understood simply as ordinary semiconductors with unusu- 
ally small gaps and strongly correlated electronic quasiparticles. An alternate interpretation 
of the resistivity, susceptibility, specific heat, and neutron scattering data is that at high 
temperatures (T » E g ) these materials most closely resemble a heavy fermion metal above 
its coherence temperature, and as the temperature drops a pseudogap opens in the density of 
states at roughly the coherence temperature as estimated from high temperature properties. 

A natural explanation of the existence of the narrow-gap semiconductors is that the 
high temperature band structure is semiconducting, and as the temperature is lowered, the 
effective hybridization of the f- levels with the conduction band is reduced as strong electronic 
correlations develop, just as in heavy electron metals. The reduced effective hybridization 
then leads directly to a dramatic reduction of the gap from its bare value. This picture is 
supported by slave boson calculations for the Anderson lattice Hamiltonian with infinite on- 
site repulsion U, for T below the bose condensation temperature ||. The mean field approach 
of those calculations does not include the finite quasiparticle lifetime resulting from strong 
electron-electron scattering. We have previously reported calculations of the self-energy 
and magnetic susceptibility for the symmetric Anderson lattice model in 2D with finite U, 
including the quasiparticle lifetime [|J. From the susceptibility and DOS we concluded that, 
on account of the temperature dependent quasiparticle lifetime, a pseudogap formed with 
decreasing temperature at a characteristic temperature much smaller than the hybridization 
gap in the bare band structure. Jarrell et al. reported essentially exact calculations for the 
D = oo symmetric Anderson lattice model with finite U ||. They also observe a metallic 
high temperature state and the formation of a pseudogap with decreasing temperature. 

Here we present additional results and a detailed analysis of the formation of the pseu- 
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dogap in the symmetric Anderson lattice model with U ~ W/2, where W is the bare 
conduction bandwidth; our results in 2D and 3D are qualitatively similar. The Anderson 
lattice Hamiltonian is ||, 

H = ~t J2 (°tcja + c) a C ia ) + J2 { V (ft c iv + Ciaficr) 
<i,j>a icr 

+ e f n{ a - h a(nl + n{ a ) } + U £ , (1) 

i 

where Cj CT annihilates a conduction (f) electron at site i with spin a, t is the nearest- 
neighbor conduction electron hopping energy, ef is the f- level energy, V is the onsite hy- 
bridization energy between conduction and f electrons, U is the Coulomb energy of two 
f-electrons on the same site, and h is a homogeneous magnetic field coupled to the electron 
spin. For the symmetric model \x = 0, and e/ is chosen so that there are exactly two elec- 
trons per lattice site at all temperatures. The chemical potential is then at the middle of 
the indirect hybridization gap in the bare band structure, between the T and M points. The 
bare activation energy Ao is half the hybridization gap, and given by 

A = ^{tDf + V 2 - Dt . (2) 

The fully renormalized single particle dispersion relations determined from the poles of 
the retarded Green's function show a temperature dependent gap E g (T), which at high 
temperatures is close to the bare hybridization gap. The high-temperature density of states 
does not have a gap, however, because the quasiparticle width at the band edges, T(T), is 
large compared to the temperature dependent activation energy A(T) = E g (T)/2. In this 
sense the system is effectively metallic (though with a poor conductivity). As the tempera- 
ture is lowered, increasing correlations in the metallic state renormalize A(T) and the gap 
decreases, but T(T) decreases even more rapidly than A(T) until T eventually reaches a 
crossover temperature T x where T(T X ) ~ A(T X ). At this point the pseudogap begins to open 
in the single particle density of states. The opening of the pseudogap feeds back into the 
calculation of the self-energy and leads to a more rapid decrease of T(T) and a rapid growth 
of A(T); these in turn accelerate the reduction of the density of states in the pseudogap 
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with decreasing temperature. At the lowest temperatures, A(T) ~ Ao- The opening of the 
pseudogap is reflected in the temperature dependence of the spin susceptibility, which drops 
faster than is possible for a semiconductor with a fixed gap when T < T x . We also note 
that for T > T x the self-energy is to an excellent approximation independent of momentum, 
but as T is reduced through T x , the self-energy becomes increasingly momentum dependent. 
This is reflected in a strong flattening of the quasiparticle dispersion at low energy. With 
decreasing U, T x increases and the explicit temperature dependence of the opening of the 
pseudogap becomes more important; the temperature dependences of the thermodynamic 
and single particle properties are less dramatic and begin to resemble those of ordinary 
semiconductors. A crossover temperature T x can still be identified from temperature depen- 
dence of the self-energy, and insofar as T X (U) is meaningful for all U, we expect that for 
T X (U) » Aq, A(T) will not differ much from A and the system will behave essentially 
like an ordinary insulator. Here we focus on the more interesting cases T X (U) < Aq. 

We have carried out fully self-consistent calculations for the 2D and 3D symmetric Ander- 
son lattice Hamiltonian using a self-consistent perturbation theory known as the fluctuation 
exchange approximation (FEA) JT2J. Here we present a highly abbreviated description of 
the FEA as applied to the Anderson lattice model and refer the reader to Ref. J7| for a 
more complete description. In a self-consistent perturbation theory, the fully renormalized 
temperature Green's function G is related to the self-energy E through Dyson's equation, 

G{k,e n ) = [G^(k,e n ) - E(k,e n )] _1 - (3) 
The self-energy is obtained from a generating functional $[G], 

and the grand thermodynamic potential is given by || 

Q(T, fi,h) = - Tr [EG + ln(-G 1 + £)] + , (5) 

which the appealing property of being stationary with respect to independent variations in 
G and E. Here 'Tr' denotes a generalized trace over all arguments of the Green's function 



and Go is the Green's function of the non-interacting system. With the spin-quantization 
axis along the field, the Green's functions and self-energies are spin diagonal, and because 
the Hubbard interaction acts only between f-electrons, $ is a functional of the f-electron 
Green's function Gf alone and only the f-electron self-energy is nonzero. This greatly 



simplifies the solution of Dyson's equation for G [IT]. The practical application of self- 
consistent perturbation theory entails the selection of a specific (possibly infinite) subset 
of the diagrams for $. The diagrams that define the FEA are those that generate the 
Hartree-Fock and second-order self-energies together with contributions from the exchange 
of longitudinal spin and density fluctuations, transverse spin fluctuations, and singlet pair 
fluctuations. Explicit expressions for these contributions to $[G] for the Anderson lattice 
model appear in Ref. 0. For the results presented below, the contribution to $ from spin 
fluctuations is within 10% of that from the second order diagram; contributions from pair 
fluctuations and from density fluctuations are smaller by factors of roughly —1/5 and —1/10 
respectively. 

We calculate Ey using an algorithm originally developed for the Hubbard model [|13|| , on 
a 32 D lattice with spot-checks on lattices as large as 128 x 128 in 2D, and using a frequency 
cutoff at least six times the unhybridized conduction band width, W = ADt. 

The magnetic susceptibility Xs{T) was found from the slope of Tr[<7 2 G] versus h for small 
fields (typically 0.005t and O.Olt); this calculation of Xs{T) is equivalent to solving an integral 
equation for the response function with irreducible vertices comprising particle-particle and 
particle- hole bubble-chains and Aslamazov-Larkin-type diagrams ||14|| . 

The retarded self-energy on the real axis was obtained from N-point Pade approximants 
|15|1 . This self-energy was then used to construct the spectral functions and quasiparticle 



dispersion relations along symmetry directions in the Brillouin zone. The densities of states 
were found by analytic continuation of the k-summed imaginary-frequency Green's functions 
using a varient of the N-point Pade approximant method of Vidberg and Serene |L6] . These 
densities of states satisfy sum rules to better than 1% accuracy. At some temperatures in 
the transition region between the metallic and insulating states (see below), it is difficult 



to obtain reliable analytic continuations. For example, Fig. 1 compares three different Pade 
continuations for the DOS at two temperatures: T = 0.055t, where the gap is opening, and 
0.035t where the gap is well formed. The three analytic continuations for T = 0.035t are 
essentially indistinguishable for all frequencies. At low energy, the analytic continuations for 
T = 0.055t agree well except very near e = 0; for It < e < 8t, two of the three continuations 
nearly coincide while the third shows qualitatively different behavior. To investigate this 
variability, we computed the rms deviation of the original imaginary frequency data from 
the data generated by using the candidate analytic continuation as a spectral function. For 
this purpose we use only frequencies less than half the high frequency cutoff to avoid cutoff- 
dependent artifacts. For T = 0.035t, the curves shown have rms deviations of 1.3 x 1CT 4 , 
1.1 x 1CT 4 and 1.3 x 1CT 4 for iV = 57, 59, and 61 respectively. At the temperature where the 
pseudogap is opening, qualitative differences in the structure of the DOS at energies ~ W/2 
are often weakly reflected in the rms deviation as seen for the T = 0.055t DOS shown in 
Fig. 1. There the rms deviations are 5.4 x 10~ 5 , 5.4 x 1(T 5 and 6.1 x 10" 5 for N = 35, 37, 
and 39 respectively; the qualitatively different continuation shows only a slightly larger rms 
deviation. The DOS results presented below will include only those temperatures at which 
the analytic continuations for a given (reasonable) N agree well with those over a range of 
N and have the smallest rms deviations. In all cases the rms deviation is less than 2 x 10~ 4 . 

For the calculations in 2D reported here, we used /i — 0, V — t and if = 0, where 
if — ef + \U(nt + n{) is the (h = 0) Hartree-renormalized f-level. With this choice of 
parameters the bare activation energy Ao is 0.2361 For the 3D model, we took V = 1.5t 
so that the ratio of the bare hybridization gap to the conduction electron bandwidth is the 
same in 3D as in 2D. We found that in this case the spin susceptibility and density of states 
for the symmetric Anderson lattice model in 3D may be roughly scaled to those in 2D by 
scaling all energies with respect to the bandwidth. An example of this scaling for the spin 
susceptibilities Xs for U = W/2 is shown in Fig. 2. In the following we present primarily 2D 
data. 

As the temperature is decreased through T x , a pseudogap opens in the single particle 
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excitation spectrum. Pseudogap formation is reflected in Xs and in the density n(/i, T); 
analytic continuation is not needed to calculate either of these. Fig. 3a shows Xs{T) for 
various U including U — 0. The form of Xs{T) for finite U strongly suggests the existence of 
a gap that apparently decreases with increasing U. For the largest U shown, U — At — W/2, 
the rapid decrease in Xs(T) occurs at a temperature ~ A /4. For U > 0, Xs(T) cannot be 
scaled to the U = result by a simple renormalization of the gap. The dependence of n 
on fi also shows a clear signature of pseudogap formation, as shown in Fig. 3b. For U = 
and T > Ao, dn/dfi is constant over — 0.3t ^ fi ^ 0.3t whereas for T = 0.0625t, dn/dfi 
shows substantial curvature. As the temperature is reduced, the gap in the single particle 
excitation spectrum becomes more apparent and dn/dfi essentially vanishes for ji < A . In 
contrast, for U = W/2 and T = 0.0625t, dn/d/i is essentially constant for —0.3 < fi < 0.3, 
consistent with a metallic state or a semiconducting state with a gap ^ 0.0625t. At the lower 
temperature of T = 0.055t, n(/i) shows slight curvature. As the temperature is lowered, 
dn/dfi — > for a region of /i about fi = as expected in the presence of a pseudogap. At 
the lowest temperature, the region where dn/dfi nearly vanishes is about the size of the bare 
hybridization gap. From n(fi) and Xs{T) one estimates essentially the same characteristic 
temperature T x for pseudogap formation. 

The formation of the pseudogap can be seen explicitly in the single-particle (tunneling) 
density of states shown for various temperatures in Fig. 4. For T = 0.125t and T = 0.03125t, 
Fig. 4a shows the total DOS over a wide energy range, from which one sees that the density 
of states is smoothly varying over most of the bandwidth. For T = 0.125t there is no 
sign of a gap at low energy, but for T = 0.03125 a pseudogap centered at e = is clearly 
evident. Over the narrower range —t<e<t shown in Fig. 4b there is still no evidence of 
a hybridization gap at the highest temperatures, but as the temperature is lowered, a sharp 
asymmetric peak signals the pile up of states at the pseudogap edge. The density of states 
at zero energy decreases rapidly with decreasing temperature. We have shown previously 
that this decrease closely tracks that of Xs(T) P|,|T7|]. 

To illustrate how the pseudogap forms, we plot in Fig. 5a the spectral functions as 
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a function of e at the T and M points which lie on either side of the pseudogap. At high 
temperature, the spectral functions are very broad due to strong electron-electron scattering. 
The tails of the two spectral functions, centered close to the bare gap edges, extend well 
into the gap where they overlap and form the observed maximum. As the temperature is 
decreased, the spectral functions at the T and M points sharpen as shown in Fig. 5b. In 
contrast to the relatively strong T dependence of the width of the spectral function T(T), 
the position of the quasiparticle band at the T point A(T) (as determined from the peak 
of the spectral function), also shown in Fig. 5b, is weakly T dependent until T(T) ~ A(T), 
which sets the crossover temperature T x . At this point the dramatic rapid decrease in 
Xs{T) and the low energy DOS begins, with a correspondingly rapid renormalization of the 
quasiparticle energy at the gap edge. 

The correlations above T x are apparent in the large renormalization of the band structure 
at low energy; the peak positions of the spectral functions at the gap edges at T x are at 
about A /4. The sharpening of the spectral functions with decreasing temperature above 
T x reflects the evolution toward coherent quasiparticle excitations. The quasiparticle renor- 
malization factor provides another measure of the strength of correlations. In a metal, the 
quantity Im E(k, e ) /e = V k {T) can be used to construct a reasonable estimate of the quasi- 
particle renormalization factor a k = (1 — <9Re£(k, e)/de\ £= o)~ 1 ', for a point on the Fermi 
surface of a Fermi liquid, V kF {T) tends to <9Re ^(kp, e)/de\ £=0 as T — > fl8f . For the one 
band Hubbard model at 1/4 filling, we found that T> kp {T) is weakly temperature dependent 
and approaches <9Re E(k, e)/de\ £=0 [fL3| . For the asymmetric Anderson lattice model 0, we 
found that the magnitude of T> k (T) increases strongly with decreasing temperature (as did 
our measure of the quasiparticle effective mass) down to the lowest temperatures that we 
studied. Here we observe that for T > T x the temperature dependence of T> k (T) resembles 
that of T> kp (T) for the asymmetric model. We plot T> k (T) for the symmetric model and for 
k at the T point in Fig. 6, in which a distinctive feature at T = T x is apparent. For T > T x , 
the system is metallic and the magnitude of T*y{T) increases with decreasing temperature. 
As the temperature is further decreased, the opening of the pseudogap acts to reduce con- 



tributions from low energy excitations to the fluctuation propagators, leading to a sharp 
reduction in the magnitude of T>-p(T), which is reflected in the widening of the renormalized 
gap shown in Fig. 5b. In the limit T — > 0, Im S(k, ie n ) tends to a finite value at low T. This 
is not inconsistent with the formation of the pseudogap; the U = Green's function already 
contains information on the bare hybridization gap, so E(k, e ) need not show a divergence 
as observed, e.g. in the Hubbard model [|19[ . 

Effects of strong correlations are also apparent in the quasiparticle band structure de- 
termined from the poles of the single particle Green's function, as shown in Fig. 7. At 
sufficiently high temperatures, the quasiparticle band structure does not differ significantly 
from the bare band structure and S(k, e) is to a good approximation independent of k. As 
the temperature is lowered toward T x , a pronounced flattening of the bands around the T 
and M points occurs and the self-energy becomes increasingly k dependent. In Fig. 7b we 
compare the band structure obtained by taking S(k, e) = S(r, e) for all k with that obtained 
from the full k-dependent self-energy. The momentum dependence of the self energy clearly 
provides a substantial contribution to the flattening of the band at low energy. Mean field 
calculations in the limit U — > oo find a quasiparticle band structure given by the expression 
for the bare band structure with a renormalized f-level and hybridization. In Fig. 7b we also 
show the band structure for a noninteracting system along the T — X direction for e/ = and 
V = 0.6t, which reproduces the particle-hole symmetry and gap of the fully renormalized 
band structure. It is clear that this approach is unable to describe the band flattening and 
band width along this direction in the FEA. 

We have presented the results of a fully self-consistent calculation of the self-energy and 
spin susceptibility in the fluctuation exchange approximation for the symmetric Anderson 
lattice model. From a metallic state at high temperature, a pseudogap forms with decreasing 
temperature. The pseudogap begins to open at the temperature T x at which the renormalized 
band gap and the width of the spectral function at the gap edge are comparable. As T is 
decreased through T x , the density of states at low energy decreases rapidly. The opening 
of the pseudogap is observed in results obtained directly from imaginary frequency axis 



data and through analytic continuation to the real frequency axis. The large temperature- 
dependent renormalizations of the band structure above T x are similar to those we have 
observed in the asymmetric model and suggest that the evolution of a coherent heavy fermion 
metallic state with decreasing T is cut off by the loss of low energy excitations as the 
pseudogap develops. 
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FIGURES 

FIG. 1. Comparison of three analytic continuations using N-point Pade approximants for 
T = 0.0552 and 0.0352 (a) low energies, < e < 2; and (b) high energies, 22 < e < 82. For 
T = 0.0552 (0.0352), 35, 37 and 39 (57, 59 and 61) Pade coefficients were used. Note that the "high 
energy" T = 0.0352 results were shifted by 0.05 for clarity and that 35 coefficients were used for 
the dashed T = 0.0552 curve. 

FIG. 2. The uniform static spin susceptibility \s for 2D and U = At = W/2 (•), 3D and 
U = 62 = W/2 (A), and the 3D result with all energies scaled by the ratio of the bare conduction 
bandwidth in 2D to that in 3D bar(o). 

FIG. 3. (a) The spin susceptibility Xs(T) for U = 0, (dashed) 2(D), 2t(o), 32(+) and 4t(«). (b) 
The total density n as a function of n for U = At and T = 0.0624t(A), 0.055t(+), 0.03125t(o), and 
0.015625(o). Also shown for comparison are n(/i) for U = and T = 0.24t (solid), 0.0625t (long 
dashed), and 0.015625£ (short dashed). 

FIG. 4. (a) The density of states over a wide energy range for T = 0.125t (short-dashed) 
in the metallic state and for T = 0.035£ (long-dashed) where a developed pseudogap is evident, 
(b) The density of states at low energy for temperatures T = 0.015625£ (long-dashed showing 
gap), 0.035t (short-dashed showing gap), 0.0625t (long-dashed-short-dashed) , 0.1252 (dashed), and 
0.25 (dotted) showing the formation of a pseudogap from a high temperature metallic state with 
decreasing T. 

FIG. 5. (a) The spectral function at the V point (solid) for T = 0.252, 0.1252, 0.06252, and 
0.0352 in order of decreasing peak width. Also shown is the spectral weight for the M point (dashed) 
and the sum of the M and V point spectral weights (long dash) for the temperature T = 0.06252 
at which the pseudogap begins to open for U = At. (b) The width of the spectral function at half 
maximum on the high energy side A and low energy side o as a function of T for k at the zone 
center, and the position of the peak of the spectral function (•) as a function of temperature. 
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FIG. 6. The function V k (T) (see text) for k at the T point and for U = t (o), 2t (□), 3t (o), 
and 4t (•). Note the sharp structure marking the opening of the pseudogap and the shift of this 
structure to lower temperature with increasing U . 

FIG. 7. (a) Fully renormalized quasiparticle band structure along high symmetry directions in 
the zone for U = 4t with T = 0.0625t (solid) and T = 0.25 (dashed) compared with the bare band 
structure (long dash-short dash), (b) Fully renormalized band structure for T = 0.0625£ (solid) 
compared with the band structure assuming a momentum independent self energy S(k, e) = S(r, e) 
(dashed), the bare band structure (short-dashed) and a fit to the mean field band structure with 
a renormalized V = 0M and ej = (long dash-short dash). 
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